arXiv: 1501.04485vl [physics.plasm-ph] 19 Jan 2015 


Linear wave action decay entailed by Landau damping in 
inhomogeneous, nonstationary and anisotropic plasma. 

Didier Benist 

CEA, DAM, DIF F-91297 Arpajon, France. 

(Dated: January 20, 2015) 

Abstract 

This paper addresses the linear propagation of an electron wave in a collisionless, inhomogeneous, 
nonstationary and anisotropic plasma. The plasma is characterized by its distribution function, 
fu, at zero order in the wave amplitude. This distribution function, from which are derived all 
the hydrodynamical quantities, may be chosen arbitrarily, provided that it solves Vlasov equation. 
Then, from the linearized version of the electrons equation of motion, and from Gauss law, is derived 
an envelope equation for the wave amplitude, assumed to evolve over time and space scales much 
larger than the oscillation periods of the wave. The envelope equation may be cast into an equation 
for the the wave action, derived from Whitham’s variational principle, that demonstrates the action 
decay due to Landau damping. Moreover, the Landau damping rate is derived at first order in the 
variations of the wave number and frequency. As briefly discussed, this paper generalizes numerous 
previous works on the subject, provides a theoretical basis for heuristic arguments regarding the 
action decay, and also addresses the propagation of an externally driven wave. 
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I. INTRODUCTION 


The linear propagation of an electrostatic wave in an inhomogeneous and nonstationaty 
plasma is, clearly, an old problem in plasma physics, that has led to numerous theoreti¬ 


cal studies for several decades (see Ref. [lj and references therein for a rather exhaustive 
bibliography on the subject). Recently, it has attracted some renewed interest, due to its 


fl. 


it is shown 


relevance to issues related to plasma compression. In particular, in Ref. 
that Whitham’s variational principle Q does apply to describe the linear propagation of 
an electron wave in an inhomogenous and nonstationary plasma. This result is actually 
not straightforward because, when the plasma density, n, varies in space, or in time, the 
Lagrangian density introduced in Whitham’s theory depends on three fields, the wave am¬ 
plitude, E p , its eikonal, tp, and the plasma density, n, while the Lagrange equations are 
only valid for E p and p>. Therefore, it is not obvious that the variations in n would entail 
variations in E p as predicted by Whitham’s theory, i.e., such that the wave action, 

be conserved. In Eq. (ED), oj = —d t <p is the wave frequency, while y is the electron suscepti¬ 
bility, defined precisely by Eq. (1641) of Section UTTl In the fluid limit, dx/dcu scales as n -1//2 


Q. 


This 


so that action conservation entails that E p should scale as n 3 / 4 , as explained in Ref. 
scaling was confirmed numerically in Ref. { 3 ], where it was shown that plasma compression 
did induce an increase in the wave amplitude, so as to keep the wave action constant. How¬ 


ever, it was also shown in Ref. Q] that compression makes the wave phase velocity decrease 
compared to the thermal one, which eventually entails the decay of the wave amplitude 
(and also of the wave action) because of Landau damping. Nevertheless, to the best of 
the author’s knowledge, there is no unambiguous theoretical proof of the latter result. In¬ 
deed, in spite of the abundant literature on the subject, no kinetic derivation of an envelope 
equation, accounting for Landau damping, could be found for an electron wave propagating 
in an inhomogeneous and nonstationary plasma. For example, the theoretical results of 


Ref. 


0 


relic on fluid equations and, therefore, they cannot account for Landau damping. 
Moreover, the expression found for y in Eq. ((TJ) is the fluid one, which becomes inaccurate 
at large temperatures. The most complete kinetic derivation of an envelope equation we are 
aware of, for an electron wave propagating in an inhomogeneous and non stationary plasma, 
is that given in Ref. [ 4 ]. This derivation stems directly from Vlasov equation, linearized 
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about the distribu 
the results of Ref. 


don function, /#(x, v, t), at zero order in the wave amplitude. However, 
4] lack of generality because fn obeys a force-free Vlasov equation, which 


is usually not the case. Moreover, and more importantly, Landau damping is not recovered 


4]. An envelope equation that does account for collisionless dissipation (and, 


m 


in Ref. 

particular, Landau damping) has been derived in Refs, jf], 6] starting from first principles. 
The equation found in these papers is valid in the linear and nonlinear regimes, and, in the 


former one, it does show t 


the derivations of Refs 


• B 


rat the wave action decays at the Landau damping rate. However, 


6 ] only hold for a homogeneous and stationary plasma. 

The present paper aims at generalizing the aforementioned articles. The main result is 
Eq. ( 1791 ) showing that, in the linear regime, the wave action decays at the Landau damping 
rate. It is derived directly from Gauss law and from the electrons equations of motion, 
linearized about the dynamics that makes the plasma evolve on hydrodynamical scales. 
This dynamics may include any kind of force field, including a magnetic one (which was not 
the case in previous publications). Moreover, the expression of the Landau damping rate we 
provide, given by Eq. ( 1701 ) , is valid at first order in the time and space variations of the wave 
number and frequency. Hence, it does not stem from a WKB approximation but vindicates, 
a posteriori , this approach. 

Among all the possible applications of these results, the present work was mainly moti¬ 
vated by the modeling of stimulated Raman scattering (SRS). Indeed, SRS is still an issue 
for inertial fusion since large Raman reflectivities were measured at the Nation Ignition 

n fin 

Facility [7j. As discussed in several papers (e.g. 8, 9]), correctly modeling collisionless dissi¬ 
pation is crucial to accurately predict Raman reflectivity. This might be done in an effective 
way by making use of envelope equations, since they proved in Refs. jlO, 11] to provide 
estimates for SRS reflectivity as accurate as a particle in cell or a Vlasov code, within much 
smaller computation times. However, the results of Refs, lfj, 11] were only for a plasma 
with constant density. There is, currently, no kinetic modeling of SRS in an inhomogeneous 
and nonstationary plasma, and the present paper is the first of a sequence of forthcoming 
articles that aim at filling this gap. Therefore, in order to make the application to SRS 
straightforward, we actually address here the propagation of a wave that may be driven. 

The derivation and presentation of our main results are organized the following way. 
Section [Hi provides a general expression of the charge density in terms of the plasma distri¬ 
bution function, f H , at zero order in the wave amplitude. This charge density appears to be 
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the sum of two contributions of distinct nature. The first contribution varies on hydro dy¬ 
namical scales, and is related to slowly-varying force held that may exist in the plasma and 
may make the density and temperature vary. The second contribution is at the origin of the 
electron plasma wave (EPW), and is used to define what may be viewed as a generalized 
electron susceptibility, 5. In Section IIII1 A(S) is expressed in terms of the variations of the 
wave number, of the frequency, and of the electric held amplitude which, from Gauss law, 
provides the wave envelope equation in a one-dimensional geometry. This equation is hrst 
derived when there is no resonant electron in Paragraph IIII A1 In this case, it is shown to 
guarantee the wave action conservation, with y in Eq. (pQ) explicitly expressed in terms of 
f H . A wave equation that accounts for the contribution of all electrons, resonant and non 
resonant, is derived in Paragraph IIII ESI Then, the expression for y is slightly changed, and, 
more importantly, it is shown that the action decays due to Landau damping. Moreover, 
our expression for the Landau damping rate accounts, at hrst order, for the space and time 
variations of the wave number and frequency. Section HVl generalizes the results of Section Hill 
to a three-dimensional geometry, while Section [V] summarizes and concludes this work. 

II. THE ELECTRON CHARGE DENSITY 

In order to describe the propagation of an electron plasma wave that may be externally 
driven, one needs to calculate the charge density, p w , induced by the EPW and the drive. 
As a hrst step of this calculation is derived, in this Section, an expression for p w in terms of 
the electron distribution function, /#, at zero order in the wave and drive amplitudes. For 
the sake of clarity, the plasma is hrst assumed to be one-dimensional (1-D), and the results 
obtained in 1-D are then easily generalized to a three-dimensional (3-D) geometry. From p w 
is also introduced, 5, defined by Eq. (TTSli or Eq. ()3TTh that may be viewed as a generalized 
electron susceptibility, and that allows to express Gauss law in a very simple fashion. 

A. One-dimensional geometry 

Let us denote by —e the charge of an electron, by m, its mass, and by f(x,v,t ) the 
electron distribution function, such that 

r+oo 

/ f(x,v,t)dv = n e (x,t), (2) 


4 






where n e is the electron density. Let us, moreover, assume that there is only one ion specie 
(the generalization to multi-ion species is straightforward, but would require unnecessarily 
complicate notations), and let us denote by rii(x,t) the ion density and by Ze its charge. 
Then, the total plasma charge density is, 

/ +oo 

f(x,v,t)dv + Zerii(x,t). (3) 

-OO 

From Liouvillc theorem, the electron distribution function is conserved, so that, 

f(x,v,t) = f 0 [x 0 (x,v,t),v 0 (x,v,t)}, (4) 

where f 0 is the initial distribution function, and x 0 and v 0 are, respectively, the initial 
position and velocity of an electron located at x, with velocity v, at time t. Let us now 
express x and v in terms of xq and Vq the following way, 


x = x 0 + v 0 t + Sxh + $x 

= X 0 + 5x (5) 

v = Vo + <5^# + Sv 

= V 0 + Sv, (6) 


where Sx and Sv are induced by the EPW and the drive, while Sxh and Svh evolve on 
hydrodynamical scales. The two latter quantities are specifically defined by, dSn/dt = Svh 
and dSvn/dt = F H /m , where F H is the force that makes all the hydrodynamical quantities, 
such as he local density and temperature, vary. More precisely, the distribution function, 
fn(X 0 , Vo,t), at zero order in the wave and drive amplitudes, satisfies the following Vlasov 
equation, 


Then, if at time t 
into, 


dfn , T/ df H , F h dfu _ n 

T" Kq—-1-7TTT - — U. 


( 7 ) 


dt dXo m dV 0 
0 , the wave and drive amplitudes are vanishingly small, Eq. (J7J) translates 


Ih(X o, Vq, t) — fo(X 0 — v 0 t — Sxhj V 0 — Svh)- 


( 8 ) 


From Eq. (JJJ), and the definitions Eqs. (J5D and (JFJ) of X 0 and V 0 , Eq. (JSJ) is equivalent to, 


f(x,v,t ) = f H (X 0 ,V 0 ,t). 


(9) 


5 





Plugging Eq. (J9J) into Eq. (J3J) for the charge density, and making use of the change of variables 
v — y Vq in the integral of Eq. ([3]), yields, 


r »+00 


p=-e 


f H (X o, Vb, t)dV 0 + Zem(x, t) - e 


f + °° dSv 

/ fH(X 0 ,V 0 ,t)-^ 7T dV 0 . 

' —oo 


dVr 


( 10 ) 


0 


Since we aim at deriving an envelope equation for the EPW that is only at first order in the 
variations of the density, we now use fn(X 0 , Vo, t) ~ fn(x,V 0 ,t ) — 5xd x fH{x,V 0 ,t). Note 
that latter approximation is not a linearization with respect to the wave amplitude, but an 
expansion at first order in the space variations of /#. Then, 


r»+oo 


"+OO 


P ~ -e / f H (x, V 0 , t)dVo + Zerii(x, t) + e Sxd x f H (x, V 0 , t)dV 0 

J— oo J — oo 

r + °° ddv f + °° ddv 

f H {x,V 0l t)—dV 0 - e / d x f H (x,V 0 ,t)Sx—dV 0 
ovo j_ QO oVo 


~ 6 

J — OO 

= Ph + Pw, 


0 


( 11 ) 

( 12 ) 


where 

/ +oo 

f H (x,V 0 ,t)dV 0 + Zerii(x,t), (13) 

-OO 

evolves on hydro dynamical scales, and is considered as a given quantity. As for the charge 
density, p w , induced by the wave and the external drive, from now on, it is identified with its 
linear value. Hence, by making use of an integration by part for the third term of Eq. (fill) , 
one finds, 

/ +oo p+oo 

Sxd x f H (x,V 0 ,t)dV 0 + e / Svf H (x,V Q ,t)dV Q , (14) 

-OO J — OO 

where we have denoted, f' H = dy 0 fH- 

We now assume that the electric held of the EPW and the drive writes, 

E = E p e i4, + zEde^+W + c . c . ? (15) 

where c.c. stands for the complex conjugate, and where E p and Ed are, respectively, the 
real amplitudes of the EPW and drive fields, that vary much more slowly than the eikonal, 
<fi. From (f) are defined the wave number, k = d x (j) and wave frequency, oj = Both 

E p and Ed are at the origin of the charge density, p w , so that E will henceforth be termed 
the total electric held, although the force Fh may include a slowly varying electric held. 
Clearly, E, may also be written the following way, 

E = E 0 e lip + c.c., (16) 
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where E 0 = ^ E p + E^ — 2E d E p sin(<5</>) is the total field amplitude (which is clearly real), 
and e* v3 = E^e^^Ep + iE d e lS ^). We moreover assume that the wave is nearly on resonance 
with the drive, so that k ~ d x p and oj ~ —dtp. 

Similarly, we write the charge density, p w , as, 

Pw = Poe* 9 * + c.c., (17) 


and introduce, 


ipo 

e 0 kE 0 ’ 


(18) 


which may be viewed as a generalized electron susceptibility. 

In the total field E, only the EPW electric field is self-consistently related to p w , since 
the drive is imposed externally (for example, in the case of SRS, the drive is due to electro¬ 
magnetic laser fields, and is just the component of the ponderomotive force along the wave 
number). Then, Gauss law is just d x (E p e u ^ + c.c.) = p w /e o, which reads, 


(■ ikE p + d x Ep)e l<t> = —ikEE 0 e l<p 

= -ikE(Ep + iE d e iS ^)e^. (19) 


Since Eq, E p and E d are real, Eq. (1T2D is equivalent to the two following equations, 

(1 + E r )E p = E d [EiCOs(5<f)) + E r sin(50)], (20) 

k~ l d x E p = E { Ep + E d [E r cos(50) - Sjsin(50)], (21) 


where 
Ref. 


12 ], 


= 3?(S) and S* = Q'(S). We henceforth assume that E d -C E p which, as shown in 
is the case for SRS. Then, Eq. (1201) is approximated by 


1 + — 0 , 


( 22 ) 


which actually yields the EPW dispersion relation. The validity of the approximation that 
leads to Eq. (122]) is discussed in detail in Ref. [jjj], where it is shown that accounting for the 
drive amplitude in the dispersion relation is important to calculate the nonlinear frequency 
shift of the plasma wave, but not to correctly estimate the linear frequency. Since we restrict 
here to linear wave propagation, Eq. (1221) may be considered as exact. 

As will be shown in Sections UTT1 and HVl 5, is of the order of (fcL 0 ) _1 or (cnTo)^ 1 , where Lq 
and Tq are, respectively, the typical length and time scales of variation of Eq, k or oj. Hence, 


7 





for the slow variations considered in this paper, Sj < 1, so that, from Eq. ()22j) . 5* <C H r . 
This lets us conclude that Eq. (l2Tf may be approximated by, 

EiE p - k~ 1 d x E p = E d cos(5(j)). (23) 


B. Three-dimensional geometry 


The results obtained in Paragraph III Al straightforwardly generalize to a three-dimensional 
geometry. Just like for a 1-D plasma, let us introduce the distribution function, 
/ff (Xq, V 0 , t), at zero order in the wave amplitude, that obeys the following Vlasov equation, 


dfn , j df H t F h df H _ n 

~df + Yo -W 0 + ~^-W 0 ~°- 

Then, Eq. (TTOl) for the charge density becomes, 

<9(V 0 + Sv) 


(24) 


p= -e 



/ff(X 0 ,V 0 ,t) 


dVn 


dV 0 + Zerii(x, t ), 


(25) 


where |<9(V 0 + 5v)/<9V 0 | is the Jacobian of the change of variables V 0 + Sv —» V 0 . The 
linear value of this Jacobian is, 


d(V 0 + Sv) 


dV 0 


lin 


i+E 

n 


dSv r , 


dVr 


On 


(26) 


where Sv n and Vo n are, respectively, the n th component of Sv and V 0 . Therefore, following 
the lines of Paragraph III Al one straightforwardly finds that the linerarized charge density 
induced by the wave and the drive is, 


Pwfct) = e 


Jx.<9 x ///(x, V 0 , t)dV 0 + e 


Sv.&v 0 f H (x., V 0 , t)dV 0 . (27) 


Let us now write p w and the total electric held as in Paragraph 111 Al 


Pw = Po e%<p + c.c., 

E = E p e^ + E d e^ +S(p) + c.c. 


= E 0 e iip + c.c., 


where E p , E^ and E 0 are real vectors, and let us introduce, 

~ _ iPo 
“ _ £ok.E 0 ' 


(28) 

(29) 

(30) 


( 31 ) 
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Then, with the same approximations as in Paragraph III At Gauss law yields the following 
equations, 


(32) 

(33) 

III. WAVE EQUATION IN A ONE-DIMENSIONAL GEOMETRY 

A. Plasma with no resonant electron 


l + ^ - 0, 


V.E ; , 

k.E„ 


k.E d 

k.Er, 


cos (S(p). 


For the sake of clarity, we first derive Eq in a one-dimensional geometry. Moreover, we 
start by assuming that there is no resonant electron, the exact meaning of this assumption 
being clarified after Eq. (|46|) and at the end of Subsection III1B 21 

From Eq. (fT4lh we now need to calculate Sx and Sv in order to derive Eq. This is done 
by solving the electrons equations of motion, 

d5x 


= 8v, 
ax 

d8v e 

= - E^ + c.c. 


(34) 

(35) 


dt m 

Eq. (I35j) is solved by making use of an integration by parts, which yields, at first order in 
the variations of Eq, 

—m 


-5v «£ 0 (t)/i(t) - ^ 


hit) + c.c., 


t'=t 


where £o{t) = Eo[x(t),t], and 


hit) = / Soit'yvw^dt’ 


(36) 


(37) 


is one primitive of £oe lip , and 

hit') = J hit')dt' (38) 

is one primitive of h- Since we only look for a first order envelope equation, we now express 
h at first order in the variations of k and c a, and h at zero order in these variations. Then, 
one easily finds, 


Hf) 

hit') 




—i 

d\b / dt ' 


£#(*') 


idip/dt') 2 ’ 


d^/dt' 2 ' 
idh/dt') 3 ’ 


(39) 

(40) 
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where we have denoted ifj(tt) = (p[x(t'),t']. 

Let us now provide an expression for Sv that is linearized with respect to Eq, the zero- 
order motion being defined by the hydro dynamical force, F H . To do so, we make use 
of the approximations, S 0 (t') ~ E 0 [x — SX 0 (t'),t'] and -0( t') ~ ip[x — 5X 0 (t'),t'], where 
8X 0 (t') = f t , V 0 (t")dt" with dV 0 /dt = F H /m. This yields, 


where, in Eq. 
Ends, 


d£ 0 /dt'\ t > =t — d t E 0 + V()d x E 0 , 
d^\) j di!\ t ’=t = — UJ + kV 0 , 

d 2 ^/dt' 2 \ t '=t = —d t u + 2V 0 d t k + V^dlk + kF H /m, 

, we used d t k = — d x co and dV 0 /dt = F H /m. From Eqs. 


(41) 

(42) 

(43) 

one easily 


h(t) 


hit) 


0 i<p(x,t) 


kV o — u 

e iv(x,t) 


+ 


d t u — kF H /m — 2V 0 d t k/dt — V^d^kfdx 1 
{kVo — a ;) 3 


(kV 0 -uj) 2 

Plugging Eqs. ilTTli and (l45l) into the expression 


(44) 

(45) 


for Sv yields, Sv = Svoe llf + c.c., with 


m 

- 5v 0 

e 


—iEn d t E 0 + V 0 d x E 0 d t u - kF H /m - 2V 0 d t k-V$d x k 

H-77T7-~- E 0 - 77T7 -^-. (46) 


kV 0 — lu ' ( kV 0 — uj ) 2 (kV o — cu) 3 

Note that our expansions, leading to Eq. (146|) for Sv, only make sense if, whatever V 0 , 
(uj — kVo)T 0 1 and (a;/Vo — k)L 0 1, where T 0 and L 0 are, respectively, the typical time 
and length scales of variation of Eo, k or uj. The latter conditions are precisely what we 
mean by the conditions for no resonant electron in the plasma. 

Since, in the expression Eq. (TTD) for p w , Sx is multiplied by d x fu , it only needs to be 
estimated at zero-order in the variations of Eo, k and u. Then, it is easily found that 
Sx = Sx 0 e lip + c.c. with, 


m 


-E, 


o 


- SX0 « —— - 

e (kv o — to) 2 


(47) 


From the knowledge of Sv and Sx given by Eqs. (1461) and (l47lh 5 is easily derived by making 
use of the expression Eq. (fl4|i for p w and of the relation Eqs. (fI7|) and (fl8|) between p w and 
5. The corresponding calculations are detailed in the Appendix |A] where 5 is expressed in 
terms of the electron susceptibility, y, defined here as, 


r»+00 


X = 


ti. 


H 


£ 0 mk kV o — u 


dVo. 


(48) 
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It is then found that S r = y, so that Eq. (fl9|) reads 


1 + X = o. 


( 49 ) 


As for Sj, we find, 


2 El L 


d_ 

dt 


dx 


E: 


d_ 

dx 


dx 


El 


X dE 0 
kE 0 dx 


ldx 

k dx 


(50) 


dco^°J dx \dk 0 

From Eq. (l49j) . 1 + x — 0 whatever x , so that d x x = 0. Moreover, we henceforth use the 
condition Ed <C E p to replace, in the expression of Sj, E 0 by E p . This actually amounts to 
deriving an envelope equation at zero order in the variations of the drive amplitude, which 
is known to be enough to correctly address SRS |10j. Then, Eq. (1251) simply reads 


!*p 2 )=2Wo s( m. 


(51) 


When Ed = 0, Eq. (T5TT) is exactly the equation one would derive from Whitham’s variational 
principle, that guarantees the conservation of the wave action defined by Eq. (111). 


B. Plasma with all electrons 

Let us now address the general situation when the plasma is composed on two kinds of 
electrons, the resonant and the non-resonant ones. In Paragraph III A1 the non-resonant 
electrons were defined by the condition, |V) — ry,| > AVr, with A; AVr To 1 and 

(A Vn/v^kLo 1. In the remainder of this paper, we use a more specific definition. 

Indeed, all the calculations are henceforth performed by assuming that, for each local value 
of the phase velocity, v Vo for the resonant electrons exactly varies between v<f, — AVr 
and + AVr (while Vo for the non-resonant electrons lies in the complementary interval). 
However, there is a technical difficulty related to this local definition for the resonant and 
non-resonant electrons, that is discussed in the Appendix [Bl As shown in this Appendix, 
this difficulty may actually be ignored (as will be done here) without affecting the wave 
equation. 


1. Basic assumption on the distribution function 

We, henceforth, assume that the resonance width AVr is such that, AVr -C Vt, where 
vt is the typical scale of variation, in velocity, of f' H about v^. For a Maxwellian, Vt would 
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just be the electron thermal velocity. For any bell-shaped distribution function, centered 
at Vo ~ 0, it is clear that there is a significant fraction of resonant electrons only if Vt is 
sufficiently large compared to v$ — A Vr. Hence, it is quite natural to assume AVr -C Vr 
when the contribution to the charge density from the resonant electrons does matter. One 
noticeable exception to the previous rule is the situation when there is a very cold electron 
beam, with mean velocity close to v^, in the plasma. Clearly, the present paper cannot 
address that situation. However, a very detailed study of the cold beam-plasma instability, 
in the linear and non linear regimes, will be the subject of forthcoming articles. 

Therefore, we shall henceforth assume AVr <C Vr, and express 5 as the sum of two 
contributions, 

^ ^res _|_ ^nres 

where S res and S nres are, respectively, from resonant and non-resonant electrons. 



2. Contribution from the resonant electrons and Landau damping 


The charge density, p™ s , induced by the resonant electrons is given by Eq. (TT4|) . with 
fu assumed to be arbitrarily small when \V 0 — v</,\ > AVr. Because of the hypothesis, 
A Vr <C Vr, to leading order, 


i-v^—AVr 

Pw s ( x , t) ~ el Sxd x [f H (x, v#, t) + (Vo - v^f'nix, v^, t)} dV 0 

J v^—AVr 

pV^+AVr 

+e / Svf' H (x,v^t)dV 0 . 

Jv^—AVr 


(53) 


For the sake of clarity, the detailed derivation of p\ff , accounting for the slow variation of 
k and oj, is reported in the Appendix O In this Subsection, we shall restrict to a homogenous 
and stationary plasma, in order to discuss in a very simple fashion the derivation of Landau 
damping. For a homogeneous and stationary plasma, 

fV^+A Vr 

Pw S {x,t)^e j Svf H (x,V'p,t)dV 0l (54) 

J V^ — A Vr 


so that p r ^ s (x,t) = p 0 e l ( kx + c.c., with, 

g 2 rt rVj+AVR 

Po = - f' H (x,v h t) / / 

m Jo Jv^-AVr 


Eolxit^AY^-^J'-^dVodt'. (55) 


12 


We now express E 0 (x,t ) in terms of a Fourier integral, E 0 (x,t ) = E K (t)e lKX dn, to find 


t />+oo 


Po ~ - — f' H {x,Vfrt) 



rv^+A V R 


e i[(k +K )V 0 -u,](t'-t) dV i d 


(56) 


E K (t')e” 

' 0 J-oo Jv^—AVr 

We, now, make use of the change of variables, u = (fc + n)AVn(t' — f), in the time integral 
of Eq. (|5gD , to hnd, 


r*+oo <r) Akx 


Po ~- f' H ( X i V<t>, t) 

m 


. . 2e*"'" ~ r ... . . Tr .shim , , , . 

/ / --E K [t + m/(/c + k)AVr]- dndu. (57) 

'-kAV R t J-o o ^ + k m 


Now, Eo is supposed to vary very little over one wavelength, so that E K is only significant 
when k -C fc, and one may assume (A: + k)AVr ~ PAVr. Moreover, it is well know that 
sin u/u is a function whose width is of the order of 2n. Hence, the integrand in Eq. (1571) 
assumes significant values only when u/(k + k)AVr < 27t//cAVr. Since, by definition, 
1/fcA Vr is much less than the typical time of variation, T 0 , of Eq, we conclude that E K [t + 
u/(k + k)AVr] may be replaced by E K (t ). Hence, at first order in the space variations of Eq, 

f + °° E K (t) 


^+oo ^ikx 


k + k 


E K [t + u/(k + k)AVr]cIk 


' — OO 
/»+00 


k + k 


e lKX dn 


K 


1 - - ) E K (t)e lKX dK 


Eq(x, t) i dE 0 


Hence, from Eq. ()57l) . we hnd, 
2 e 2 


Po 


mk 


f H (x,v^,t)e^ x -^ 


E 0 (x,t) + 


k 


i dEn 


k 2 dx 


k dx 


Sill u 


du. 


'-kAV It t u 


(58) 


(59) 


We shall now restrict to times, t, such that kAVRt 1, a condition that will be discussed 
in great detail in a few lines. Then, in the integral of Eq. (l5Tj) . — kAVRt may be replaced by 
—oo, so that 


Po 


7re 

mk 


f'llfaVfrt) 


Eq{x, t) + 


i DEq 


k dx 


(60) 


From Eq. (1601) and the definition Eq. (ITT)) of 5, we conclude that, for a homogenous and 
stationary plasma, 


27TC 


-J' H {x,V^,t) 


1 + 


i dE 0 
kE 0 dx 


(61) 


Eomk 2 ‘ 

As is quite obvious, and will be shown in Subsection IIII B 11 this term, combined with the 
contribution from the non-resonant electrons, will provide the Landau damping of the wave 
action. 


13 























Let us now discuss physically the condition kAVut 1 (in practice kAVut > 2 tt is 
enough), we had to impose in order to approximate po by the expression given in Eq. (j60]h 
As we shall now show it, this condition is actually not restrictive. Indeed, we only discuss 
the propagation of an EPW once its amplitude is much larger than the noise level (otherwise, 
it makes non sense to assume slow variations for E 0 ). Therefore, this EPW has necessarily 
grown, either by itself or due to some external drive, at least during a finite time. Hence, 
if we denote by t g the typical time of variation of the EPW during its growth, we only 
consider times t t g . For a wave that keeps growing, t g is the same as To, the typical 
timescale of variation of E 0 . In this case, the condition kAVut 1 is necessarily fulfilled 
since, by definition, kAV^T 0 1. However, one may also consider the situation when, 
after having been driven for a finite time, the EPW experiences Landau damping. Then, 
during the damping phase, T 0 ~ 1/z/r, where z/r is the Landau damping rate. In this case, 
AVr is such that kAV^/vi 1, and Eq. (1601) for p 0 is on ly accurate after a time t such 
that kAVRt > 2n. Hence, our expression for S res , and therefore our derivation of Landau 
damping, is only valid after a time that scales as 1 /z/r, although it may be small compared 
to 1 /z/r. Now, it is well known that there exists a finite time before the settlement of Landau 
damping, because Landau damping is only an asymptotic result. We, therefore, conclude 
that the condition kAVut 1, which we had to impose in order to obtain the expression 
Eq. (l6Tf for S res , is not restrictive. 

Moreoever, it is quite interesting to discuss further, in terms of AVr, the conditions under 
which the value of po, given by Eq. ()57lh is close to its approximate expression, Eq. ()60lh 
There are basically two conditions. The first one is that the integral, f^ kAV//t sin u/udu, has 
to be close to ir/2, which requires a resonance width, AVr > 27 r/kt. The second one is that 
E K (t + u/kAVu) ~ E K (t), which requires a resonance width that has to be, at least, of the 
order of 1/T 0 . When T 0 ~ 1 /z/r, we recover the result of Refs. [l3, 14] that the resonance 
width for Landau damping first decreases as 1/t, until it saturates to a value of the order of 
vi/k. Here, this result is generalized to a wave that may grow or decay, and whose amplitude 
may vary in space and time. 

The initial decrease of the resonance width, AVr ~ 2n /kt, has a very simple physical 
meaning, that we shall now explain. We define AVr such that, when \V 0 — > AVr, 5v 

may be accurately approximated by Eq. (l46|h This approximation does not only require 
a nearly linear electron motion but, also, a nearly adiabatic one. Indeed, the first term in 
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Eq. (1461) is at zero order in the space and time variations of the fields and, actually, it is 
easily found that this term is only the linear limit of the adiabatic value of 5 derived in 
Ref. 15]. The next two terms of Eq. (UH may be simply viewed as first order terms in 


an expansion about a purely adiabatic motion. Hence, the non-resonant electrons are the 
nearly adiabatic ones. Now, as is well known, only after a finite time, t a , may one neglect 
the non-adiabatic part of the electron response. Here, we just find that t a ~ 27r/fcAI/R, that 
is the the period of an electron orbit in the wave frame (provided that Vo is very different 
from v^). 

Since the non-resonant electrons are the nearly adiabatic ones, we conclude from 
Section IIII Al that, when all electrons are nearly adiabatic, the wave action is conserved, 
as predicted by Whitham theory. Hence, the conservation of the wave action, defined 
by Eq. o. is closely related to the conservation of the electron dynamical action, 
/ = {k/2n) J Q 27r//fc vdx. This is not surprising since Whitham theory introduces a local 
Lagrangian density, which, therefore, assumes a local electron response, that is only true if 
the electrons are adiabatic. 


In order to conclude this Subsection, let us now report the results of the Appendix [Cl 
where S res is derived for a nonstationary and inhomogeneous plasma. In particular, account 
is taken of the space and time variations of the wave number and frequency. Then, Eq. (16T|) 
is changed into, 


i dk i dE 0 
k 2 dx k dx 




£q mk 2 


Ih (Xj V<fi 7 f). 


(62) 


3. Contribution from the non-resonant electrons 


The contribution, S nres , from the non-resonant electrons, has already been derived in 
Section IIII Al Using the results of this Section, we introduce, 

/*+oo r'nres 
JH 


X = - t / - dV 0 , (63) 

e 0 mk J^ kV 0 - uj 

where fff es is the distribution of the non-resonant electrons, fH res (x,Vo,t) & 0 when |Vo — 
Vtpl < AVr. Since AVr is assumed to be much less than the typical scale, vt, of variation of 
f Hl one may approximate x by, 


X 


£ 0 ink 


P.P. 


n 


H 


kVn — uj 


dV t n 


(64) 
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where fn now accounts for all electrons, and where P.P. stands for the Cauchy principal 
part. Now that x is defined, one may use the results of Section 1111 Al for the real part of 


= X, 


and the expression Eq. (|50|) for its imaginary part, 


2 Eg L 


d 


dt\du °) dx\dk 0 


X dE 0 


ldx 

k dx 


(65) 


( 66 ) 


kE 0 dx 

Using Eqs. (|62li and (j65h . one finds that the equation 1 + S r = 0 [i.e., Eq. (|22li of Para¬ 
graph OTX] deduced from Gauss law] reads, 


X = -i + 


7re 


£ 0 mk 2 


f H (x,v^t) 


1 dk 


+ 


1 dEn 


k 2 dx kE 0 dx 


Hence, up to terms of order two in the space variations of the fields, 


1 


2EI L 


S f d\ 


dt \du 


d (d\ 


- T- ^E 


dx \ dk 


+ 


1 dE 0 
kEn dx 


(67) 


( 68 ) 


4■ Landau damping of the wave action 


In order to find the equation ruling the evolution of the wave action, we first need to use 
the same approximation as in Paragraph III Al that, in the expression for Eq deduced from 
Eqs. (|S2I> and (EEJ); one may replace E 0 by E p . Then Eq. (123]) . ’E i E p — k~ 1 d x E p = E d cos(<5</>), 


yields, 


where 


d ( d X 9 (dx „ 2 \ O.. d X tti 2 


m\8£ e -) ~m\dk E d + 2VL d£ E -~ 2 E p e ^ os ^ 


7re J 


El = 


v h*)< 


(69) 


(70) 


£„mk 2 d w x' 

is the Landau damping rate. It is noteworthy that the expression for z/^, at first order in 


the variations of k and a?, is exactly the same as for a homogeneous and stationary plasma, 
which vindicates a WKB approach for Landau damping. 

With our definition, f H is proportional to the local electron density, n e (x,t), and one 
may want to introduce the local plasma frequency, u> pe = n e e 2 /e^m, in order to provide a 
more familiar expression for vl, 

— f (7A\ 

L ~ k 2 d u x n e (x,t ) 

Clearly, Eq. (l69j) shows that, if the wave is not driven, its action decreases at twice the 
Landau damping rate. 


16 






















IV. WAVE EQUATION IN A THREE-DIMENSIONAL GEOMETRY 


The derivation of the wave equation in a three-dimensional geometry proceeds as in 
Paragraph IIII E>[ where the contributions to the charge density from the resonant and non¬ 
resonant electrons are calculated separately. 

The expression for the charge density induced by the resonant electrons is a plain gener¬ 
alization of the 1-D result Eq. (|53j) . namely, 


Pu 


+e 


pv^—AVr 

/ <fx.<9 x V j_, t ) + (14 - v^)dv k fH kj dV k ) dV± 

'v^—AVr 
rv^+A V R 


<5v .dv 0 fH\v&y ± dVk <7V_l, 


(72) 

Iv^-AVr 

where 14 is the projection of V 0 along the local direction of the wave vector, and Vj_ is 
perpendicular to that direction. The derivation of p r ^ s is detailed in the Appendix O where 
it is found that the value of S res is exactly the same as in 1-D, namely, 


—i 


. i dk i dE k 

1 - — — + 


7re 


£o mk 2 


/h(x,i'*,<), 


(73) 


k 2 dx k k dxk 

where, x k and E k are, respectively, the local direction and the component of the field am¬ 
plitude along k, and where, now, we have denoted, f' H = ff dv k fffdV±. 

As regards the contribution from the non-resonant electrons, it is expressed in terms of 
X, which is defined as in 1-D, 


X = 


"+oo f’nres 

Jh 


£ 0 7Tlk ,/_ 00 kV k ~ UJ 

/ r f ; 


dV k 


£ 0 7Tlk 


P.P. 


H 


kV k — uj 


dV k 


Then, as shown in the Appendix |Aj the real and imaginary parts of 


(74) 


are, respectively, 


^nres 
l—J r 

^nres 
1—l z 


= X, 

d 2 x _ v g kX + k - (3kX-V)Ecil - Ep.Vx - x'V.Eq 


(75) 

(76) 


Then, Eq. 


from Gauss law, 1 + S r = 0, yields, 


, 1 dk 1 dE k 

X = -i - —— + 


(77) 


k 2 dx k ' k dx k 

so that, at first order in the space variations of the fields, Vy ~ 0 and —yV.E 0 ~ V.E 0 . 
Moreover, we make use of the same approximation as in 1-D and replace, in the expression 
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for Hj, E 0 by E p . Then, Eq. fl33jl from Gauss law, E t — V.E p /k.E p = k.EdCos(<5</>)/k.E p , 
reads, 

d Lx ~ V4x , k. [duxdtEp - (dkX- V)E P ] vre 2 , k.E^ 

-2- + --*,*,() = —cosW). (78) 

Now, for a plasma wave, up to terms which are of order one in its space variations, E p & 
Ek(k/k) (see, for example, Ref. jic|]). Hence, up to terms which are of order two in the fields 
variations, one may use in Eq. (JZHJ), E p fa Euik/k) ~ E p (k/k), which yields, 

I (§H - (I*?) + = 2E ^ cos «>- < 79 > 

where the Landau damping rate, v L , is still given by Eq. (1701) or (171]) . Hence, one recovers 
the 1-D result that, if the wave is not driven, its action decays at twice the Landau damping 
rate. 

V. CONCLUSION 

Starting from first principles, we described, in this paper, the linear propagation of an 
electron plasma wave. Our analysis accounted for any type of force held that might make the 
plasma evolve, in space and time, on hydrodynamical scales. Therefore, our main assumption 
was to consider the plasma as collisionless. However, we also had to exclude the very special 
situation when there was, in the plasma, an extremely cold beam whose mean speed was the 
wave phase velocity. Except for the two former hypotheses, our analysis was totally general, 
and came to the very simple conclusion that, if the wave was not driven, its action decayed 
at the Landau damping rate, Moreover, we found no correction to the WKB expression 
for vl, at first order in the variations of k and to. Therefore, we showed that the concept of 
wave action was very robust, valid in the most general situation. Moreover, we provided a 
solid theoretical basis to the simple use of a WKB formula in order to predict the evolution 
of the wave action, a result that was not completely obvious, although abundantly used. 

When deriving the Landau damping rate, in the envelope equation for the EPW, we 
showed that the electrons that mainly contributed to it were those whose velocity lied in the 
range, [v^ — AV,ry, + AH], where AV was time-dependent. For short times, AV ~ 1/kt, 
until it reached max(2n/kT 0 ,2™$/kL 0 ), T 0 and L 0 being, respectively, the typical time and 
length scales of variation of the EPW amplitude. This extended previous results on the 
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resonance width for Landau damping [l3, 14] to an EPW amplitude that may vary in space 
and time, and that does not necessarily decays. 

Our analysis of the EPW propagation was made general enough to account for an external 
drive, which was necessary to be able to address the important issue of stimulated Raman 
scattering, that motivated this work. Actually, the present paper is the cornerstone of 
forthcoming articles that discuss how to model SRS by accounting for kinetic effects in an 
inhomogeneous and non stationary plasma, in the linear and non linear regimes. 

As for the linear propagation of an EPW, this is, clearly, one of the oldest issue in plasma 
physics, if not the oldest one. As stated in the introduction, it was not obvious at all that, 
for any kind of plasma, one could simply model it by the blind application of Whitham’s 
theory, completed by the blind use of a WKB formula for the Landau damping rate. This 
result is now proved. 
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Appendix A: Derivation of when there is no resonant electron 


1. One-dimensional geometry 


In order to derive the electron charge density, p w , given by Eq. (fl~4lh we use the results of 
Section UlI A1 that Sv = Sv 0 e z,f + c.c. and Sx = Sx 0 e llf + c.c., with Sv 0 and Sx 0 respectively 
given by Eqs. (1461) and (l47lh Then, from Eq. (ITSli is obtained an expression for 5, in terms 
of the electron susceptibility, y, which, when there is no resonant electron, is 


X = 


f+CXD ft 

JH 


-dV 0 . 


(A.l) 


£ 0 mk J_ 00 kV o — co 

Now, from the expression Eq. (1471) of Paragraph III Al for Sx, and by making use of an 
integration by parts, one finds, 

r*+oo 


/»+oo ^2 

e Sxd x f H {x, V 0 ,t)dV 0 = 


e z E 0 e^ f + °° d x f' H 


■dV o + c.c. 


(A.2) 


_ mk .Woo kVo — co 

Using Eq. flA.2(l . together with the expressions Eq. (|46|) for Sv, Eq. (TT4]) for p w yields, 
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Pw = Po£ ltp + c.c, with 
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(«Vb - w) 

We now take advantage of the fact that /# obeys the Vlasov equation © to find, 
(. F H /m)f' H = —dtfn — Vod x fn- This yields, by making use of an integration by parts, 


r»+00 


k(F H /m)f' h 
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Therefore, from Eqs. (1A.3D and (1A.4I) . one easily finds that 5 = ip 0 /e 0 kE 0 is, 
i d\ dE 0 


(A.4) 
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(A.5) 


£q mk J _ OQ kV o — uj 
where we made use of the identity dtk = — d x uo. 

Now, it is clear that the second line of the right-hand side of Eq. (1A.5I) is just 
(i/2)d t (dx/dou), while the third and fourth lines of the right-hand side of Eq. (IA.5[) are, 
respectively, —(i/2)d x (dx/dk) and —( i/k)d x x ■ This straightforwardly yields, 


= X, 


2 E 2 L 


dt\du °) dx\dk 0 


J^dEo 
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1 dx 

k dx ’ 


(A.6) 

(A.7) 


where 5, = 'A(S) and 5 r = 3?(S). 


2. Three-dimensional geometry 


In a three-dimensional geometry, we find it more convenient to make use of an integration 
by parts in the second term of Eq. (1271) to find 


Pr(xd) = e 



d5v r, 


<Sx A/S™(X, V„,i)dV 0 -e JJJJ2 gvpj»“^ V ». «)<Wo, (A.8) 
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where Sv n and Vo n are, respectively, the n th component of 5v and Vo- Similarly, by making 
use of an integration by parts, we find, 


X = 


£o m 



fn 


(kV 0 -uj) 


M. 


(A.9) 


Recalling that the total electric field is, E = E 0 e llf + c.c., and using the same method as 
in Paragraph IIII A1 one easily finds that, to the relevant orders, 5x = fixoe l<f + c.c. and 
<5v = dv 0 e lip + c.c. with, 
eE 0 


dx 0 = 
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Now, even if the plasma is magnetized, the n th component of F# does not depend on the 
n th component of Vq. Therefore, 
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By making use of an integration by parts along 14 (the component of V 0 along k), and 
taking advantage of the fact that k.F h does not depend on 14, one finds 
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where Fj_ and Vi are, respectively, the components of F# and V 0 transverse to k. Then, 
since Fj_ does not depend on Vj_, from Eqs. (1A.13I) and f lA.141) . 
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Moreover, taking advantage of the identity <9 t k = —Vu, one finds, 
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Therefore, plugging Eqs. (lA.IOj) and (1A.11I) into the definition Eq. (1A.8|) for p^ res , and using 
the identities Eqs. flA.15j) and (1A.16H . one finds that p™ res reads, p™ es = p 0 e lip + c.c., with 


po _ k.E 0 
£o 2 

Hence, the real and imaginary parts of 


[dlx - V4x] +k. [d uX dtE 0 - (a kX .V)Eo]-E 0 .Vx-xV.Eo-^k.E 0 . (A.17) 
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Appendix B: Definition of the resonant and non-resonant electrons 

In Sectionini the non-resonant electrons are defined to be such that | Vo— v^\ > AVr. This 
definition needs to be global, i.e., for the non-resonant electrons, Vo G (—oo, V\) U (V 2 , + 00 ), 
where V\ and V 2 are fixed, independent of space and time, and such that, whatever the local 
value of Vfp, — V\ > AVr and V 2 — > AVr. Indeed, if we defined the non-resonant 

electrons such that Vo G (— 00 , — AVr) U (v# + AVr, + 00 ) then, the boundaries of the 

integral in Eq. (1631) for % would depend on v<j,. Consequently, the derivatives of x with 
respect of k and uj would involve the derivatives of these boundaries, thus affecting the 
results derived in Paragraph IIII Al 

Using a global definition for the non-resonant electrons entails that, for each value of 
Vff,, the velocity Vo for the resonant electrons is such that Vo G [v^ — AV\ , + AV 2 ], where 

min(AVi, AV 2 ) > AVr. Then, from Eq. (l56lh the contribution, p r ^ s , to the charge density 
from the resonant electrons writes p^ es = pg es e lip + c.c., with 
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The first term in Eq. (1C.ID . involving the velocity integral from v$ — V to + V , clearly 
yields the value for pg es derived in Paragraph IIII B 2l since AVi > AVr. As for the second 
term, involving the velocity integral from + V\ to v^, + V 2 , it is proportional to 
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with AV m = (AVi + Ai/ 2)/2 and AV = (AVi — AV 2 V 2 . Eq. (jC.2j) shows that J 2 oscillates 
with time at a period close to T = 2ti/ kAV m . Therefore, the value for p r u f s derived in 
Paragraph IIII B 21 may be viewed as the averaged value of p r ^ s over T, so that the envelope 
equation is also an averaged equation over T. ffowever, since AV m > AVr, and since 
AVr is chosen so that kAVR T 0 , where T 0 is the typical time of variation of Eq, k or oj, 
an envelope equation averaged over T = 2n/kAV m may be considered as the instantaneous 
one. 

Note, moreover, that if AV AV m , J 2 is of the order of Eq x m.&x.(AV/AV m , T 0 //cA14n)- 
In this case, the correction to p r u f s is negligible at any time, and not only on the average. 


Now, it may happen that the phase velocity is so different in two remote locations that 
the conditions \Vq — v,p\ > AVr over the whole plasma, and AVr <C vt, nray not be both 
fulhlled. In this case, we simply divide the plasma into domains of finite length. These are 
such that the conditions |Vo — v<f,\ > AVr over the whole domain, and AVr vt may be 
both fulhlled. Then, for each domain, the envelope equation (1691) is valid, with x very close 
to the value given by Eq. 064[) . Hence, the same envelope equation holds over each domain, 
so that this equation is valid for the whole plasma. 


Appendix C: Derivation of H res for a nonstationary and inhomogeneous plasma 
1. One-dimensional geometry 

In 1-D, the value of E res is derived from that of the charge density, p™ s , induced by the 
resonant electrons, and given by Eq. (153|) of Paragraph (IIII B 2(1 . 

p r w es = Pi + P2, (C.l) 

where 

rv^—AVii 
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rv#+AV r 

P 2 ~ e Svf' H (x,v < p,t)dV 0 . (C.3) 

J v^—AVr 

Let us start by estimating the term p x . Since we restrict to hrst order equations in space and 
time, we only need to express 5x at zero order in the variations of Eq, k and 00 . Moreover, 
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the force Fh may be neglected in the electrons dynamics, so that, using the notations 
£ 0 (t) = E 0 [x(t),t], one finds, 


8x « — e i ( kx - ut) £ Q {t")e i{ - kVo ~ u:)(t "- t) dt"dt' + c.c. 


t rt' 


m 



(C.4) 
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Therefore, using the fact that 1/kAVn is much less than the typical time of variation of £q, 
one finds that, when A V^t 1, 


rv</>—AV r 

/ Sxdv o ~ - 

S>-A V R m 


2 F e i' 





—me 


o J o 
t rt' 


k(t — t") 


m 


e i{kx-ut) / / So{f)5{t - t")dt'di" + c.c. 


o Jo 


= 0. 


(C.5) 


Similarly, 
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= 0. 
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We, therefore, conclude that, at first order in the fields variations, p\ = 0. 


Let us now estimate p 2 , using the same method as in Subsection IIII B 2l Hence, we start 
by expressing E 0 [x(t),t)] as a Fourier integral, E 0 (x,t ) = E K (t)e tKX dK. Then, since the 
total field is, E(x,t ) = E 0 e l,f( - X F _|_ c introducing the notation 

a K (x, t) = kx + <p(x, t), (C.7) 


one finds, 


+oo 


E(x, t ) = / E K (t)e ia V x F dK + c.c. 


(C.8) 


J — OO 

Since we restrict here to the linear response with respect to E, this lets us express the 
electron position the following way, 


x(t') « x(t) + V 0 (t' 

= x(t) + V 0 (t' 


t) + 
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dt'" dt' 1 
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t) + Sx H , 


(C.9) 
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where Fr is the force term in the Vlasov equation for /#, Eq. ©, accounting for the fact 
that the plasma is inhomogeneous and nonstationary. 

Hence, from Eqs. (1C.3(1 . flQ.8j) and (1C.9ft . one finds 

—e 2 f +oa U 

p 2 = — —f' H {x,v < j ) ,t)j / E K (t')I 3 (t' - t)dt'dK + c.c., (C.10) 

^ J —oo Jo 


where 

/•v^+AVji 

I 3 (t' - t) = / e ia4x+V 0 (t'-t)+5x H ,t'] dVo ' (C.ii) 

J v$~A Vr 

The integral I 3 (t' — t ) is calculated exactly the same way as the integral I\ (t) of Para¬ 
graph HTX] which yields, 
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da K /dV 0 ( da K /dV 0 ) 3 J J v ^_ AVr 

l[x + (v# + AV r )( t' -t) + 6x H , t ']e iaK[x+{v * + A ' / «)h'-d+fe«T] 

(C.12) 
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-X[x + {vf - A V R )(t' -t) + Sx H A^ [{v ^ VR){t '~ t)+dXH,t ' ] ) , 


where 
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k(x, t) + K 
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(C.13) 


k((,t’) + K [ [k((,t') + k] 2 _ ' 

By hypothesis, AVr is much smaller than the typical scale of variation in velocity, vt, of f' H . 
Unless for very peculiar distribution functions, this entails AVr <V. v Moreover, we also 
assume that Sxr -C v^t — t'), which would usually be the case since Sx R is due to the space 
and time variations of the density, which changes on scales much larger than the fields E 0 , 
k and oo. Therefore, since I is a slowly varying function, one may use the approximation, 
X[x + (v ( j,AAVR)(t , — t) + SxH,t'] ~ I[x Fv^t' — t),t']. Moreover, we henceforth also neglect 
Sxr in the phase a K . Since this phase varies quickly, its local value does depend on Sx r , 
however, because (v<j> ± A Vr){ 1J — t ) varies much more rapidly with t' than 8xr, the time 
integral of I 3 (t' — t ) should not be much affected by the neglect of 5xr- Hence, we come to 
the following approximate value for J 3 , 


J 3 (+> _ ~ JJ[x-\-V^jt t),t] (ia K [ x+ ( Vil>+ AV R )(t'-t),t'] _ e ia K [(v < / > -AVR){t'-t),t')'\ _ (C.14) 

[k(x, t ) + n](t' — t) \ / 

Clearly, when t' —> t, I 3 [t' — t) —> AVr, while, when 1 1' — t\ > 2n/(k + k)AVr, I 3 {t' — t ) 
oscillates with t' , the amplitude of the oscillations decreasing as 1 /(k + n){t' — t ). Then, 
since E K only assumes significant values when k <C fc, and since 1 /kAV R is much less than 
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the typical time of variations of E K , Is(t' — t ) may be considered as a delta function in the 
time integral of Eq. (1C.101) . Hence, just as in Section UlIB 21 in this time integral, E K (t r ) 
may be replaced by E K (t). Then, Eq. (1C. 101) translates into, 

” 2 /,+ °° E K (t)e iKX 
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7re 

rnk 


f H (x,v h ty^ 


1 + n/k 
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where 
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, h (x,t) = — / llx+v4t'-t),t'] e -‘^ + ^- -—-f- dt'. 

71 Jo \t — t) 
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Since k -C k, we may now approximate Eq. (1C. 151) by its expansion about k — 0, i.e., make 
use of the approximations, 1/(1 + n/k) ~ 1 — n/k and rj K ~ rj + hid K r] K | K=0 , where we have 
denoted rj = r/ re=0 • This yields, 
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i dE n 


k dx 


+ C.C., 


(C.17) 


where the term d x E 0 x d K rj K | re=0 has been neglected because it is of order two in the fields 
variations. As for r] = r] K= 0 , using Eqs. (1C.7|) . (1C. 13j) and (1C.16|) . it reads, 


V = Vi + V2, 


(C.18) 
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Clearly, at first order in the variations of k and oj, 

i dk 
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In order to estimate rji, we introduce, 

5tp±(x,lf) = ip[x + ± A V R )(t' - t),t'} - ip(x,t) - k^ ± A V R ){t’ - t)] - u(t' - t) 

= (p[x + (Vtf, ± AV R )(t' ^ t),t'] - (p(x,t) T kAV R (t' - t). (C.22) 


Then, for r small enough, 
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Using the hypothesis v</, AVr, one hnds that, when r ~ 1/kAVn , 

0</> ± A V r ) 2 t 2 dk dk 


dx 


k 2 AV 2 dx 
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1 dk 


(C.24) 


kAVR k 2 dx 

Now, AVr is chosen so that uj/AVr be much larger that the typical length of variation, L 0 , 
of k. Hence, 
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Similarly, 
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(kL 0 ) 2 co 2 dt 

Therefore, over a time of the order of 1/kAVR, 5ip± varies in a negligible fashion. This 
shows that e l <p[®+(¥>±AVR)(t , -t),i'] oscillates at a frequency very close to kAVR , so that the 
integral in the right-hand side of Eq. (1C.19|) does not vary much after t > 27t//cAVr. In 
this integral, it is therefore licit to make use of the approximation 8ip + ~ ~ Sip = 

ip[x + v^it' — t),t'] — ip(x,t). Actually, from Eqs. (1C.25jl - flC.27|) . it is clear that the Taylor 
expansion Eq. (1C.2311 is valid up to a time r much larger than kL 0 /kAVR = Lq/ AVr. Since 
AVr, over such a large time, it is likely that resonant electrons have travelled over the 
whole plasma length. Hence, the same Taylor expansion may be used for 5ip + and Sip- so 
that, indeed, 5(p + ~ 8ip_ ~ Sip. We, therefore, obtain the following approximate value for 
hi, 
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Let us now introduce S(t') = (2/7r) f t , 2sin[kAVR(t" — t)]/{t" — t)dt" . Then, by making use 
of an integration by parts, one hnds, 


r/i ~ 1 — k(x, t) 


/•t e i5ip 
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Now, it is clear that S(t') is a Heaviside-like function that varies from 0 to 1 over a time 
interval of the order of 1/kAVR. Therefore, by making use of a Taylor expansion for 5p>, one 
finds, 

3(771 - 1)-r 2 [vld x k - v$ x i 0 - d t u\ , (C.30) 

where r is the order of 1/kAVR. Since, in this paper, we only account for terms which are, 
at least, of the order of k~ 2 d x k, ( uk)~ l d x io , or u!~ 2 d t uj : from Eqs. (1C.25|) - (1C.27D . we find 
3(?/i — 1) ~ 0. 

Similarly, 

9) , Ft h -\- f).h 

(C.31) 
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k 2 AV R 

Therefore, from Eqs. flC.171) . (1C. 181) and (1C.211) one hnds that, at first order in the fields 
variations, 
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k 2 dx k dx 

Using Eq. (1C.II) for p r u f s and the fact that p\ ~ 0, one hnds that p™ s ~ p 2 is given by 
Eq. (1C.32I) . Then, from the very dehnition Eq. (TlST) of S, 
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2. Three-dimensional geometry 


The charge density induced by the resonant electrons in a three-dimensional geometry is 
given by Eq. ([72]) of Section [TV] and writes 


Pw S = Pl + P2, 


(C.34) 
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’v^-AVr 

Just like in 1-D, pi is calculated at zero order in the variations of k and c 0 . Then, the 
component Sx n of Jx is given by Eq. (IC.4j) with £ 0 replaced by £ n , the component of the 
held amplitude along the n th direction, and Vo replaced by V k . Then, the integration over 
14 straightforwardly yields p\ ~ 0. 
































As regards P 2 , in 1-D we found that, except for the term (1 /k 1 2 3 4 5 6 7 )(dk/dx) in Z, the variations 
of k and c j are negligible for the following reasons, i) The velocity integral yields a delta¬ 
like function whose temporal width is of the order of l/fcAV#. ii) Over this time interval, 
the resonant electrons move by about v^/kAVn, which is much less than the typical length 
of variation, L 0 , of k and oj. Similarly, in 3-D, the integration over 14 will yield a term 
proportional to ( l/k 2 )(dk/ dx k ), where x k is along the local direction of k. Except for this 
term, the longitudinal variations of k and oj are negligible for the same reasons as in 1- 
D. As for their transverse variations, they are also negligible provided that V± < (kLo)v^. 
If we assume that, for such large values of V± (usually larger than the speed of light), 
/# (x, v#, Vx) ~ 0, then <5v may be considered as a function of k.V 0 — oj. Consequently, by 
integration over Vx, the contribution to P2 from the components of hv perpendicular to k 
yield a negligible contribution. As for the component of Sv along k, its contribution to p 2 
is exactly the same as in 1-D. 

Therefore, we conclude that the 3-D expression of S res is a mere generalization of 
Eq. (1C. 3 31) obtained in 1-D, namely, 
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where, E k is the projection of the held amplitude along k, and where, now, we have denoted, 


/ h ( x , Vfrt) 
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